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A random graph order is a partial order achieved by independently sprinkling relations on a vertex 
set (each with probability p) and adding relations to satisfy the requirement of transitivity. A post 
is an element in a partially ordered set which is related to every other element. Alon et al. [2] 
proved a result for the average number of posts among the elements {1, 2, . . . , n} in a random 
graph order on Z. We refine this result by providing an expression for the average number of 
posts in a random graph order on {1, 2, . . . , n}, thereby quantifying the edge effects associated 
with the elements Z\{1, 2, . . . , n}. Specifically, we prove that the expected number of posts in a 
random graph order of size n is asymptotically linear in n with a positive ^/-intercept. The error 
associated with this approximation decreases monotonically and rapidly in n, permitting accurate 
computation of the expected number of posts for any n and p. We also prove, as a lemma, a bound 
on the difference between the Euler function and its partial products that may be of interest in its 
own right. 



1. Introduction 

Several definitions of random partial orders can be found in the combinatorics literature. If the 
number n of elements of the underlying set is fixed, perhaps the most natural definition is that of 
"uniform random order," in which we pick a member of the set of n-element posets uniformly 
at random. Although no practical way is known of generating posets according to this definition 
for large n, it is known that as n ^ cx) most of them are "3-level" posets f6l. In this case, 
increasing n leads to posets with a greater width but not a greater height, on average, because 
of the growth in the number of relations per element. A second definition is that of "random 
fc-dimensional order", for some integer k, in which one picks k linear orders on n elements 
uniformly at random (in other words, k randomly chosen permutations of the set {1,2,..., n}), 
and then takes their intersection. The resulting posets are of dimension k by construction, and 



some of their properties are known 1 13 1. 



A third definition, and the one we will mainly be interested in here, is that of "random graph 
order" which depends on a parameter p e (0, 1). To obtain a partial order of this type, one 
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first generates a random graph on the vertex set {l,2,...,n} by including an edge with 
probabiHty p for each pair of vertices i and j; one then turns the graph into a directed one by 
converting each edge (i, j) into a relation z -< j in the partial order \f i < j (in the usual order 
on the integers); and, finally, one imposes transitive closure by adding relations so that i < k 
whenever there exists a j such that i < j and j -< k. Several properties of random graph orders, 
such as width, height, and dimension have been studied |[T] . In particular, it is known |[T] that the 
expected height of a random graph order on n elements grows linearly with n. 

In the physics literature, random graph orders are also known as "transitive percolation" be- 
cause they arise in a special case of the theory of directed percolation |7 1, where non-local bonds 
in a 1 -dimensional lattice are turned on with probability p. They also play a prominent role 
among the stochastic sequential growth models that have been proposed for the classical version 
of the dynamics of causal sets |9 |, and this is the context that most directly motivates our work. 
A causal set |4| is a partially ordered set that is locally finite, meaning that the interval or Alexan- 
drov set j) := {l \ i <l < j} is finite for every pair with i -< j. In the causal set approach to 
quantum gravity (for a recent review, see Ref. |5 1), the poset is seen as a discrete spacetime. The 
partial order corresponds to the causal relations among its elements, and "i -< f can be read as 
"i causally precedes j", while volumes of spacetime regions correspond to the cardinality of the 
appropriate subsets of the poset. The final theory is expected to assign a spacetime volume of the 
order of £p to each element, where £p := yjGhf(? = 1.6 x 10^"^^ cm is the Planck length. 

A post in a poset is an element that is related to every other element in the poset. In other 
words, each post n divides the ordered set into the subset of elements that precede it, its "past", 
and the subset of elements that follow it, its "future". In the causal set interpretation the space- 
time "pinches off at n; this can be seen as the zero-spatial-volume singularity at the end of a 
collapsing phase for the universe and the beginning of a new expanding phase. Thus, a first set 
of desirable conditions for transitive percolation to be considered as a physically reasonable way 
of generating discrete spacetimes is that if a random graph order develops multiple posts, the 
number of elements between two posts be allowed to grow sufficiently large for that region to be 
able to model our observed universe. 

It has been known for some time fZ) that infinite random graph orders have an infinite number 
of posts. However, the occurrence of posts in finite random graph orders has not been studied as 
extensively. We begin by revisiting random graph orders on N and then analyzing finite posets. 
Roughly speaking, we find in our analysis that "edge effects" are small but non-negligible. In the 
infinite case, the probability that element n is a post approaches a constant value fairly rapidly. 
In a finite case, the expected number posts is well approximated by this limiting probabiHty 
times the size of the set plus a small positive offset. We illustrate our conclusions with numerical 
simulations. 



2. Infinite Random Grapli Orders 

We begin by calculating the probability that any given element in an infinite random graph order 
is a post. In order to express this probability succinctly we define \ — p. 
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Figure 1. Plot of Afc(g) as a function of g for fc = 1, 8 (from top to bottom). 



The function hi{q) is known as the Euler function, and it has been studied in considerable detail 
p2) . In particular, K{q) > for all < (7 < 1. We now prove a theorem for one-way infinite 
random graph orders similar to the result of Alon, et al. |j^ that the probability that an element 
in a random graph order on Z is a post is ^^{q). 



Theorem 2.1. In a random graph order on N with probability < p < 1, the probability that 
any given element k is a post is given by 



PT{kisapost) = Afe_i(q) Aoo(g) 



> K^f- (2.1) 



where by >, we mean "greater than and, for large k, approximately equal to." 



Plots of \k{q) as functions of q for various values of k (see Fig. [T]) illustrate the rate of the 
convergence of Xk{q) to K{q)- Because of this convergence, the similarity in (2.1 1 holds for 
"most" k E N. This observation suggests two intuitively plausible results. The first is that with 
unit probability, there are infinitely many posts in any random graph order. This result was first 
proved by Alon et al. Q. Although the original result was for random graph orders on Z, the 
proof is easily adapted to partial orders on N. The second result is that the mean spacing between 
posts is Equivalently, the expected number of posts after N stages is well-approximated 

by k'^N. This information on the structure of random graph orders is of the type we are interested 
in from the point of view of their possible applications as models of discrete spacetimes, and in 
the next section we will consider it in more detail. 



Proof of theorem. If k is related to A; — 1, fc — 2, . . ., fc — i + 1 (for 1 < i < k), then the 
probability that k ~ i -/( fc is q*, since by transitivity the only way for this to happen is for k ~ i 
to be unrelated to each of fc — i + 1, fc — z + 2, k. Hence we find that 

Pr(fc - i^k\k-l^k A k-2^k A... A k - i + I ^ k) ^ I ~ q\, 

where we have used the notation 'Pr{A\B) for the conditional probability of A given B and A 
for logical and. Repeatedly decomposing the probability that k is related to each element before 
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it yields the expression 

Pr(fc related to all previous m) = Pr(A; - 1 -< fc) • Pr(fc - 2-</c|fc-l^fc)-... 

• Pr(l ^k\2^k A... A k - I ^ k) 

= Afc_i(9) 

> >^{q)- (2.2) 

The inequality follows because the partial products of Aqo (<z) are strictly decreasing. 

On the other hand, the same logic that led to ^2.2} shows the probability that n is related to 
every later element is given by 

Pr(?i related to every later m) = Pr(n -< n + 1) • Pr(ri -<rt + 2|n^n+l)-... 

= nr=i(i-9^) 

K(q). 

Moreover, the events m ^ n for m < n and n ^ m for m > n are independent. If transitive 
closure were to relate m < n and fc > n in a manner involving n, then n would be the middle 
element and both relations m ^ n and n ^ k would exist a priori. Hence the event "n is a post" 
will occur if and only if the two preceding, independent events occur, which has probability 

Pr(n is a post) = Pr(n related to all m < n) ■ Pr{n related to all m > n) 

= A„_i(g) Xoc{q) 

as desired. □ 



3. Posts in Finite Random Graph Orders 

From the results for infinite graph orders we expect that, to a good approximation, the expected 
number of posts in an n-element poset is -/Vposts('T-) = ^^(9) n. In fact, the mean number of 
posts among the elements {1, 2, . . . , n} in a random graph order on Z is equal to K?{q) n 
However, this number should be smaller than the expected number of posts in a random graph 
order on {1,2,..., — and appreciably smaller for small n — ^because elements near the edge 
are significantly likely to be related to all the elements in {1, 2, . . . , n} but not all the elements 
in Z. We have carried out numerical simulations of transitive percolation with different values 
of p and n; Fig. [2] shows the resulting values of A^posts plotted versus n for a fixed p. From this 
plot, one may see that for large n the number of posts is well approximated by a line with a small 
offset. We will now prove the following theorem. 

Theorem 3.1. For all < q < I, there exists a sequence of real numbers {&n(9)}^i ■so that 
for all 71 > 1, bn{q) is strictly between and 1 and the expectation value {N^^q) of the number 
of posts in a transitively percolated causal set on {1, 2, 3, . . . , n} with probability p = I — q 
satisfies 

{N„,g) = n^{q)-n + KM)- (3.1) 
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Figure 2. Average and standard deviation of the number of posts in a random graph order versus 
number of elements, for p = 0.35. For each value of n, four hundred n-element causal sets 
were generated, and the average and standard deviation of the numbers of posts are shown. The 
sampled values of n are the multiples of 5 up to 20,000. Notice in the inset graph that the average 
number of posts seems to be well fit by a line with positive y-intercept. 



Moreover, {bn{q)}'^=2 strictly monotonically decreasing to a positive limit b{q) given by the 
expression 



h{q) = 2K{q)Y,{Mq) -<<!))■ 



(3.2) 



fc=0 



For notational convenience, we will drop the explicit dependence of hP' and 6„ on q. We also 
introduce the abbreviations 



nk 



k=l 



Notice that /zi — k. In order to prove the theorem, we first establish the following estimates. 



Lemma 3.2. For all n > 2, we have 

(A„_i ~ 9")5„ < q" 

and 



(3.3) 



(3.4) 



Proof. If A„_i — g" is not positive, then (3.3 i clearly holds since 5„ and the right-hand side 



are positive. So suppose that A„_i — q" is positive. Because the Afe are monotonically decreasing 
in k and n — 1 > 1, we may replace A„_i with Ai to find 



{Xn-i-q^Sn < (Ai-g") 



^ q-^ 
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Distributing and using the definitions of Xk gives 



E 



E 



q 



n(fe+l) 



„n/c „nk 

\ ^ q \ ^ g 



- q') ■ 



where we have reindexed the second sum in the second line. Separating the first term in the first 
sum, we get 



q 



" + E 



u=2YlUi^-q') t2Yl\=A^-q') 



E 



q 



nk 



V - 



1 — g*^ \ — q 



<q" 



because all the terms in the summation are negative. This establishes (3.3 1. 

Now we will show that ( 3.4 1 follows from ( 3.3 i. We use the following well-known identity ||8) 
which holds for all complex |a;| < 1, \z\ < 1: 



11(1--"-)- E 



(3.5) 
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For completeness, we include a proof of (3.5 i. Define 

f{x,z) 



n™=i(i-^™^)' 



(3.6) 



and consider / as a function of z for fixed |a;| < 1. Since nm=i(l ^ a;™^) 7^ 0, we can write it 
as exp (XIto=i ~ x"^z)). This, in turn, may be written 



exp 



00 00 

EE 

m=l fc=l 



exp 



EE 

k=l m=l 



k 



This shows that / is analytic throughout the unit disk \z\ < 1. Therefore we can write / as a 
power series in z: f{x, z) = J2T=o Ck{x)z^ . Now it is easy to see from (3.6 1 that /(a;, xz) = 
(1 — xz) f{x, z), and this implies x'^ Ck{x) — Ck{x) — x Ck-i{x). Noticing that co{x) = 1 and 
solving this recursive relationship gives (|3.5[l. 



Finally, set x = g and z — ^, and apply ( 3.5 1 to the definition of /i„ to obtain: 

/ ^ , \ -1 



E 



Notice that ^„ = (1 + Sn) ^, and rearrange (3.3 1 to get A„_i S'„ < (1 + S'„) g". Putting these 
together, we have 

A„-l - K = A„_i(l - = , , ^ < q , 



as desired. 



□ 



Lemma 2 provides a very useful bound regarding the convergence of the partial products of 
K to their limiting value. We will use the bound several times to prove statements that involve 
expressions of the form A„_i — k or 1 — fi„. While there are more efficient ways to compute 
the Euler function [12], Lemma 2 also shows that even the naive products converge reasonably 
well: the error is strictly bounded by g", a considerable improvement — especially for q close to 
1 — over the obvious estimate that the partial products are of order 0{q^^). Slater [11 J first used 
( |3.5| l to compute the Euler function numerically but did not observe its implications for the naive 
products. 



Proof of theorem. First, we define the random variables 



1 fc is a post 
fc is not a post 



(fc = l,2,...,n). 



Then N^^g = J2k=i and, by linearity, {Nn^q) = YJl=i{Xk) = YJl=i Pr(fc is a post). From 
the proof of Theorem 2. 1 we know 

fc— 1 n—k 

Pr(A: is a post) = 

i=i 3=1 
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Substituting in the definitions of and k gives: 



fc=i 



/£ = ! 



Define the "offset" quantities 6„ according to (3.1 



k=l 



1 



- 1 



(3.7) 



First, we produce a lower bound on 6„. As x < —log (1 — a;) for a; e (0, 1), we have 



logfik = -^log(l-g*) > = 



i=k 



i—k 



which impUes that 



From this we obtain 



fj.k < exp 



1 > exp 



qk _j_ qU—k+i 



fJ-k Mn-fc+l V 1 ^ g 

Substituting the previous equation into p.7[) gives 



-1 > 



gA; _|_ gM— fc+1 



1-9 



L ^ 2 g 
Dn > K 2^ 



fc=l 



2^^2(1-9") 
(1-9)^ ' 



(3.8) 



where we have used the fact that the exponents n — k + 1 and k range over the same set of values 
as fc = 1,2, ... ,n. This establishes that the sequence {bn}^^2 bounded below by a positive 
quantity (see Fig.[3]|. 

Next, we prove that for every n > 2 we have 6„ > hn+i ■ We calculate the difference bn—bn+i : 



K - K+l) 



E 

k=l 



^J'k^J■n-k+l 



n+1 

E 

k=l 



1 



tJ'kfJ'n-k+2 



= 1 - 



1 



A*lMn+l 



" 1 



k=l 



1 



n+1 



1 



Recall that summation by parts (see, e.g., Ref |101) says that for general sequences {x„} and 
{y„}, if we define X„ = J2k=i ^k, we have ^21=1 ^kVk = X^Vn + Xk{yk - yfe+i)- 

Taking Xk = l/^in-k+i - '^/lJ-n-k+2 and yu = in this formula, we get (notice that Xk is 
a telescoping sum): 



1 / 1 



E 

A:=l 



The quantity 1 — 1/ in the first set of parentheses in the sum can be made smaller 

by replacing with /i^^ as /i„ is monotonically increasing in n. Performing the remaining 
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telescoping sum yields 

1 Wl I \ f I I \ f I 1 

> 1 
= 1 



MlMn+l \Ml P-n+l J \ fJ-n / \ Ml 

2/1 1 \ 1 



Using the fact that n — jii — A„_iMn and expressing all the denominators in terms of n gives 

= ^ [k^ - A^_i + 2A„_i - 2A„] . 
Factoring out a factor of A„_i from each variable in the numerator leads to 

= ^ [Xn-l{^in - + 1) + 2?"] . 

Now, because /i„ — 1 is negative, we can replace /i„ + 1 with 2 to make the first term in brackets 
more negative. Then we have 

K-^{hn - bn+l) > (A„_i(/i„ - 1) + g") 

= ^^(«;-A„_i +(?"), (3.9) 

which is positive by the preceding lemma. This establishes that {&ri}^2 monotonically de- 
creasing and hence must converge to a (positive) limit b. 



Now we will show that 6„ < 1 for all n. From (3.7 1, we have < 6i < 1 since hi — 1 — k^. 
Also, 62 = 2(1 — q — K^), so 62 < 1 whenever q + > 1/2. We use the inequality k > 
1- q - q^, which holds for all < q < 1 by \3A). We get that q+ > q + {l- q - q^y. 



To verify that the right-hand side is greater than 1/2 for all q, notice that its derivative factors as 
— (1 — 2(7)(1 + 4g -I- 2q^), which implies that it achieves its minimum at g = 1/2. At g = 1/2, 
it is equal to 9/16, which proves that ^2 < 1. Since is monotonically decreasing for n > 2, 
this establishes 6„ < 1 V n. 



Finally, we will prove ( 3.2 1. Define B to be the right-hand side of (|3.2[i — with Bk its fcth partial 



sum — so that our goal is to show b = B. Begin by writing H = [ri/2j and Soddin) ~ 1 if n is 



odd and if n is even. By splitting the symmetric sum in ( 3.7 1, we have 

H 



&„ = 2 ^ (Afc_iA„_fc - K^) + (5odd(?^)(A^ - k^) 

k = l 
H 



fc=i 



Now since /i^ij,_|_2^ > 1, we can replace it with 1 and take the limit as n ^ cx) (so that the 
^odd{n){Xj^ — K^) term goes to 0) to get that b > B. Notice that S is a positive series that is 
bounded above, and hence it must converge. Now we look at the difference between the partial 
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sums bn and Bh-i'- 

K - Bh-1 = Afc_i (Ai,T-fc+i ^ l) + 3odd{n){\H - k^) 



H 



fc=l 
H 



< ^ {f^n-k+i - 1) + '5odd(r^)(Al^ - K^) 

k=l 

< 2kH - 1) + 6odd{n){X% - K^); 



in the last step we used to the facts that n — [n/2\ > [n/2\ and that fik is monotonically 
increasing in k. Multiplying and dividing the first term by Xh and using the fact that Ai/(1 — 
Mff+i) < q^^^ by Lemma 3.2 yields: 

bn - Bh-1 < 2kH-^ + S^ddin)iXjj - k^) 

^2Hq"+' + 5odd{n){X%~n^). 
Taking the limit as n ^ oo of both sides gives that b < B, which completes the proof. □ 



4. Conclusions 

In this paper we considered random graph orders on n elements. Consistent with the known fact 
that infinite random graph orders have infinitely many posts [|2j, we have shown in Theorem 2 
that the mean value of the number N^^q of posts grows Unearly in n with a mean separation 
between posts approaching 

r OO -1—2 

L fe=i ^ 

the same mean spacing between posts as in a random graph order over Z. 

In a finite random graph order over {1,2,..., n}, however, the actual value of {Nn.q) is not 
exactly proportional to n but includes a small positive offset bn{q). This offset stems from the 
fact that the first and the last few elements are appreciably more likely than the other ones to 
be posts. Thus, the functions 5„(g) (shown in Fig. [sjl quantify the edge effects, with ^ bn{q) 
corresponding to the contribution of each of the two ends of the poset. As Theorem 2 and Fig. 
|3]show, that contribution does not vanish even in the n — > cx) limit. Although we do not have 
an analytical expression or bound for the standard deviation of Nn.q at this point, our numerical 
results suggest that — consistent with intuition — it is proportional to ^/n (see Fig.|2]i. 

Overall, this result does not significantly affect the number of posts or the size of the inter- 
post region in a random graph order. As mentioned in the Introduction, the latter is one of the 
first quantities one considers when estimating how viable such posets are as discrete models 
for spacetime. Therefore, as far as allowing inter-post regions of a random graph order to grow 
large, the only condition we need to impose is that q be sufficiently close to 1; equivalently, 
the probability p = 1 — q of linking two elements in the random graph must be small enough. 
The size of the inter-post regions is not the only condition we would impose for a poset to be 
manifoldlike; we are also studying the effects of other requirements and will report our results 
on them separately. One interesting byproduct of the work reported here, however, is the bound 
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(3.4 1 in Lemma 2. This bound was useful for us in proving the assertions in Theorem 3.1 but 



because of the importance of the Euler function, it may be useful in other contexts as well. 
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